function Gauss_quad_1D_trial_test(coe_fun,?)
int_value = 0;
 
for k=1:Gpn
    int_value=int_value + Gauss_weights(k)*feval(coe_fun,Gauss_nodes(k))*FE_basis_local_fun_1D(?)*FE_basis_local_fun_1D(?);
end